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Abstract. The gap probability generating function has as its coefficients the 
probability of an interval containing exactly k eigenvalues. For scaled random 
matrices with orthogonal symmetry, and the interval at the hard or soft spec- 
trum edge, the gap probability generating functions have the special property 
that they can be evaluated in terms of Painleve transcendents. The deriva- 
tion of these results makes use of formulas for the same generating function 
in certain scaled, superimposed ensembles expressed in terms of its correlation 
functions. It is shown that by a judicious choice of the superimposed en- 
sembles, the scaled limit necessary to derive these formulas can be rigorously 
justified by a straight forward analysis. 



1. Introduction 

1.1. An applied setting for gap probabilities. The first use of random 
matrices to problems in theoretical physics was in relation to the study of the 
spectra of heavy nuclei (see |Por65| for a collection of early works on the subject). 
It was hypothesized that the highly excited energy levels of heavy nuclei would 
have the same statistical properties as the the eigenvalues from an ensemble of large 
random real symmetric matrices. More explicitly, the large random real symmetric 
matrices were chosen from the Gaussian orthogonal ensemble (GOE) in which each 
matrix X occurs with probability density e~ TlX / 2 (such an ensemble is invariant 
under the transformation X i— ► OXO T where O is a real orthogonal matrix; this 
has the physical interpretation of there being no preferential basis and explains 
too the adjective orthogonal in GOE). To leading order matrices from the GOE 
have the support of their eigenvalues in [—y2N, y/2N]. The largest eigenvalue thus 
occurs in the neighbourhood of \/2N, which is referred to as the soft edge, while 
the region away from the edges (for example in the neighbourhood of the origin) 
is referred to as the bulk. It is the statistical properties of the eigenvalues of large 
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GOE matrices in the bulk, scaled so that their mean spacing is unity, which are 
compared against data from the spectra of heavy nuclei (with the latter also scaled 
so that the mean spacing between consecutive levels is unity) . 

More recently problems from statistical physics have led to applications of 
distributions at the soft edge of random matrix ensembles (i.e. in the neighbourhood 
of the largest eigenvalue.) Here each of the three symmetry classes — orthogonal, 
unitary and symplectic — are relevant to the applications. In terms of Gaussian 
matrices, such symmetry classes are realized by the probability density e~^ TtX / 2 in 
which the Hcrmitian matrix X has real elements in the case of orthogonal symmetry 
(ft = 1), complex elements in the case of unitary symmetry (/? = 2) and real 
quaternion elements represented as 2 x 2 matrices in the case of symplectic symmetry 
(/? = 4). In the latter case the corresponding matrix X has doubly degenerate 
eigenvalues, and the convention is to count only one of the distinct eigenvalues in 
the trace. For each of the three symmetry classes the corresponding eigenvalue 
probability density function (p.d.f.) is of the form 

1 N 

(i-i) ciis/sfo) n \**-*if 

1=1 l<j<k<N 

with g[x) — e~P x2 (C denotes the normalization). Independent of (3 the leading 
order support is on [— \2N, y/2N], and in the neighbourhood of the largest eigen- 
value at x = \^2N (the soft edge) the eigenvalues have spacing 0(l/iV 1 / 6 ). With 
EN,p{k; (s, oo);gp(x)) denoting the probability that exactly k eigenvalues are in the 
interval (s,oo) of the ensemble as specified by (Sjy,/? is broadly referred to as 
a gap probabilty), in keeping with these facts one expects 

(1.2) lim E N Jk](V2N + cpX/N 1/6 ,oo);e-l 3x2/2 ) 

to be a well defined order one quantity (here cp is an N independent factor chosen 
for convenience). Significantly, it is also expected that for all gp(x) such that the 
eigenvalue support is to leading order a single interval with right endpoint a(N), 
there will be a scale b(N) such that 

(1.3) lim E N Jk;(a{N) + b(N)X,oo);g (x)) 

A 7 "— too V / 

exists and is equal to (|1.2[) . Such universality questions have been a major theme 
of P. Dcift and collaborators (see |Dei99j for a summary of this work up to 1999 
relating to (3 = 2, and [DG041 IDG061 ICDG061 IDGKV06j for recent results on 
/? = 1 and 4). In particular, with 

9f3 (x) = x a e~^l\ (x > 0) 

which corresponds to the so called Laguerre ensemble, it is known from the rigorous 
work of Johansson and Johnstone |Joh00l IJohOlj that for j3 = 1 and 2 

lim E N ^(k; (V2N + X/V2N 1/6 , oo); e _/3ai ") 

= lim E Nif3 (k;(4N + 2(27V) 1 / 3 , oo); a^e"^ 2 ). 

This limiting probability is denoted E s ^ olt (k; (s, oo)). 

One class of problems in statistical physics giving rise to these probabilities 
relates to a last passage percolation problem, originally formulated by Hammersley, 
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and various symmerizations due to Baik and Rains [BRQl]. Here a unit square, 
with bottom left corner at the origin, contains points uniformly at random with a 
Poisson rate of intensity z 1 (this means that for 8 small, each non-overlapping S x S 
subsquare has probability z 2 5 2 of containing a point). Continuous piecewise linear 
paths are formed from (0, 0) to (1, 1) by joining points in the square with segments 
of positive slope. The length of a path is defined as the number of points it passes 
through, and l u is used to denote the maximum of the length of all possible paths. 
The limit theorem of Baik, Deift and Johansson |BDJ99| tells us that 

lim Prf^^ < y) = E s 2 oft (0; (y, oo)), 

so relating to the soft edge gap probability with (3 = 2. In regards to this probability 
with (3 = 1, let I s /2 denote the maximum length of all paths going from (0,0) to 
this diagonal. A limit theorem of Baik and Rains [BROlj gives 

MPt(^ < y) = ^(0; (y,oc)). 

Furthermore, with the points constrained to be symmetrical about the lower left to 
upper right diagonal, Baik and Rains [BROlj proved an analogous limit theorem 
relating to i?| oft . 



1.2. The gap probability generating function and Fredholm determi- 
nants. A special feature of the gap probabilities {Ep oit (k; (s, oo))}fc = o,i.... for the 
random matrix couplings /3 = 1,2 and 4 is that they can be expressed in terms 
of Fredholm determinants and Painleve transcendents. How this comes about has 
different features for each of the three fi values. However a common step is that 
one introduces the generating function 

oo 

En,p(J; g(x); 0=^(1- O k E N , (k; J; g (x)) 

k=0 

(because E^^{k; J;g (xj) — for k > N the sum terminantes at k — N). About 
£ = this has the expansion 

(1.4) EN,p(J;g(x);Q = 1 + „, J dx x --- J dx n P( n )(xi, . . . ,x n ) 



n=l 



where pr n \ denotes the n-point correlation function for the point process specified 
by tm 

The case f3 = 2 is the simplest, because the p.d.f. then specifies a determi- 
nantal point process, which means that its n-point correlation is an n x n determi- 
nant. Explicitly, with {Pj(^)}j=o,i,... denoting the monic orthogonal polynomials 
with respect to the weight function g2(x), and (a, b) 2 ■= f_ g2(x)a(x)b(x) dx, one 
has that (see e.g. [Forj l 

(1.5) P( n )(xi, ■ ■ ■ ,x n ) = det K N (xi,Xj) 
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where 



(1.6) 



JV-l 

K N (x,y) = (g 2 (x)g 2 (y)) 



1/2 Pj{x)Pj(y) 



(g2(x)g 2 {y)) 1/2 Pn(x)pn~i(v) - pn-i(x)pn(v) 



{PN-1,PN-1)2 x-y 

with the final equality due to the Christoffcl-Darboux summation formula. Sub- 
stituting (|1.5p in (jl.4p one recognises the final expression as an expansion of a 
Fredholm determinant (see |WW65] ). giving 

(1.7) E Nj2 (J;g 2 (x);$ = det(I- (K N>J ) 

where Km,,j denotes the integral operator on J with kernel ()1.6|) . With g 2 (x) cor- 
responding to the Gaussian or Laguerre weights as introduced above, the soft edge 
scaling limit is easy to perform rigorously JohOO , leading to the result [Fbr93j 

(1-8) E%*((s, oo);0 = det(I - ^% } ) 

where K^° { ^ is the integral operator on (s, oo) with kernel 

(1.9) *■*(*,„) = Ai( g )Ai^ y )-Ai( g )Ai-(x) 

x-y 

This is the sought Fredholm determinant evaluation, and it in turn can be used 
|TW94|, IAvM95l IBD02] to deduce a Painleve transcendent evaluation. 

The situation at (3 = 1 and 4 is more complex. The correlations are now 
quaternion determinants rather than scalar determinants. In relation to the former, 
let X be a 2N x 2N antisymmetric matrix, and set 

-1 



Z 2 N '■= IjV 



1 



With the quaternion dual of a 2N x 2N matrix specified by 

A° = Z 2 nA t 

one notes that matrices of the form XZ 2 m are self quaternion dual. The quaternion 
determinant qdet of such matrices has the key property |Dys72| that 

qdetXZ 2N = Pf X 

where Pf denotes the Pfaffian. As a consequence, for A self quaternion dual, the 
quaternion determinant and ordinary determinant are related by 

(1.10) (qdet A) 2 = detA 

Analogous to (|1.7| . one has that if, as is the case at (3 = 1 and 4 (see e.g. [ForJ), 



p( n )(xi, . . . ,x n ) = det K N (x ly Xj) 



where Kn is a 2 x 2 matrix such that [Kn{xi, £j)]jj=i,...,n is self quaternion dual, 
then 

EnA J ; 9pix); =qdet(I-^jv,j). 
Here -fC/v,j is a 2 x 2 matrix integral operator on the interval J, and i"Gv,j has kernel 
KN{x,y). Squaring both sides and applying (|1.10[) shows 

(1.11) (E N}0 {J;gp{x);O) 2 = det(I- 
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The formula with £ = 1 is taken as the starting point of the analysis 

in [TW96j leading to Painleve evaluations of £f oft (0; (s,oo)) and £| oft (0; (s,oo)). 
This strategy is used in [Die05| to give the analogous formulas for the corresponding 
generating functions. However the details of these calculations are very technical. 
Fortunately there is another approach to the problem introduced by the present 
author in ForQO], and to be further developed herein. 

In this approach, instead of working with the soft edge scaled limit of (jl.lljl , 
and thereby involving a matrix integral operator, the first step is to derive the 
alternative expression 

(1.12) (£:r ft (0;(s,oo)) =Ef t (0;(s,oc))(l-J s (I - )~ 1 A s [y]B s (y) dyj 

where K^° 1 ^ is as in (|1.8p . A s is the operator which multiplies by Ai(a;), while B s is 

the integral operator with kernel J °° Ai(y — v) dv. Starting from (|1.12| the Painleve 
expression can be deduced in a page or two of working (assuming knowledge of 
some identities from [TW94] ) . Another significant feature of (|1.12|) is that it has 
been taken as the starting point of the proof of the Fredholm determinant formula 
|Sas05L IFS05] 

Ef t (0;(s,oo) = det(l-V^ ) ) 

where V solt (x, y) is the integral operator on (0, oo) with kernel V solt (x, y) = A\(x + 
y + s). 

To extend the approach of ForOOj to the Painleve evaluation of the generating 
function B| offc ((s, oo); £) |For06] . a key identity is the formula 
(1.13) 

/•OO 

S odd(OBsoft)= ((s!0o);0 = dct(I-^°t ) )(l-^ J [(I~^Klf oo) )-'A s }(y)B s (y)dy) 

which for £ = 1 is equivalent to (|1.12[) . Here, with the ensemble (|1.1[) for (3 = 1 
referred to as OEjy(gi(x)) (here the "O" denotes the underlying orthogonal symme- 
try), OEjv(gi(x))UOEjv(gi(x)) denoting the superposition of two independent such 
ensembles, and the operation "odd" referring to observing only each odd labelled 
eigenvalue in the superposition as counted from the right most eigenvalue (i.e. the 
soft edge), odd(OEsoft) 2 refers to the soft edge scaling limit of the ensemble 

(1.14) 0dd(0E A ,(e- ;r2 )U0E w (e- ;E2 )). 

Another advantage of the approach of [ForOO is that the Painleve evaluations 
for (3 = 4 are deduced as a corollary of those at j3 = 1 and (3 = 2. This is possible 
because of inter-relations between the gap probabilities for the three symmetry 
types jFROlj . 

The primary motivating factor in seeking an alternative approach to the deriva- 
tion of the Painleve evaluations of the (3=1 and 4 soft edge gap probabilities was 
to calculate analogous formulas for the hard edge. The hard edge refers to the 
neighbourhood of the origin when gp(x) ~ x a as x — ► + , gp{x) = for x < 0. In 
the case gp{x) = x a e~ l3x / 2 , x > 0, the appropriate hard edge scaling is x i— > X/AN, 
and working based on superimposed ensembles can be carried through, leading to 
the sought Painleve evaluations |ForOO|. IFor06] . More explicitly, let 

(1.15) 0dd(0Ejv(x (a - 1)/2 e-^ 2 ) U OEjv^^^^e-^ 2 )) =: even^OEjv) 2 
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denote the joint distribution of all odd labelled eigenvalues in the superposition 
(labelled from the hard edge at x = 0). Let odd(OEhard) 2 refer to the hard edge 
scaling of this joint distribution. The identity which plays the role of (|1.13|) at the 
hard edge is 



=(o-l)/2 



£odd(OEhard) 2 ((0;S) .£. a) 

(1.16) = dct(I - C^)) (l - £ / ~[(I - e^)- x A h ] (y)B h (y) dy) . 



Here the argument a refers to the hard edge singularity x a , K^ r ^ is the integral 
operator on (0, s) with kernel 

hard. \ _ JaWx)y/yJ'aWv) ~ V^afV^Wv^) 



(1.17) K^ {x ,y) 2{x _ y) 

while A h (x,y) is the operator which multiplies by J a (y/x) and B h is the integral 
operator on (0, s) with kernel J a (t) dt. 

1.3. Aim of the paper. In this paper we will reconsider the derivation of 
(|1.13|> and (I1.16|) . which in the unpublished work }For00j is not rigorous. We will 
begin by recalling and giving a critique of the strategy used in [ForOO] . Then we will 
proceed to present a rigorous strategy based on different ensembles than those used 
in (|1.13p and (|1.16|) to scale to odd(OEsoft) 2 and odd(OEhard) 2 respectively. These 
ensembles have the advantage that the closed form expressions for the correlation 
functions are simpler than those of the original ensembles. This makes the analysis 
of the scaling of the corresponding gap probabilities much simpler. 



2. Gap probability generating function for superimposed ensembles 

2.1. Review and critique of the original calculation. To study the deriva- 
tion of (|1.13p and (|1.16p one first notes that for an integral operator I + C ® D, the 
fact that C <g> D is of rank 1 gives that 

/oo 
C(y)D(y) dy. 
-oo 

It follows that (|1.13[) and (| 1 . 16|) can be rewritten 

(2.1) £ odd ( OEsoft > 2 (( S ,oo);£) = det(l-£(i^))+i s ®S s )) 

(2.2) £°dd(OEhard) 2 ( ( 0j s );£ ( a _ i)/ 2 ) = dot (l - ^(^(ts)) +A h (g>B h )y 

Here the tilde on the operators indicates that the kernels have been multiplied by 
the gauge factor {\ / d{x))a{y) for a decaying sufficiently rapidly so that the integral 
operators inside the determinants are trace class. This latter technicality can be 
avoided by expanding (|2.1[) and (|2.2[) according to the right hand side of (|1.4p to 
obtain 

£odd(OEsoft) 2 ((S;0o) . e) 

(2.3) = 1 + j— / dxi--- dx n p°d d (° E&oft ) ( a;i) . . . 5 Xn ) 
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with 
(2.4) 

odd(OEsoft) 2 
P{n) 



(xi, . . . , x n ) = det K soit (xj,Xk) + M(xj) I Ai(xk~v)dv 



and 



£odd(OEhard) 2 ^ £. ^ _ 

" (-0" " 



n=l 



, odd(OEhard) 2 , 



(2.5) 

with 
(2.6) 

oddCOEhard)^^ . >S J = det [tfW^) + _ 



J a (t)dt 



j.k— l,...,n 



The derivation of (|2.3|) given in jForOO] took as its starting point an explicit 
form of the n-point correlations 

odd(GOEjv) 2 / \ i , \ T r i \ 

P (n) {xi,...,x n ) =det K N -i(Xj,Xk) 



(2.7) +e-^/ 2 2-( A '- 1 )i7 Ar _ 1 (^)(A 1 ( a;fe ) + A 2 {x k )) 

where K^_i is specified by (jl.6|) with ^(^) = e _a:2 , and 



j,fe=l,...,n 



, / x e_r/2 (iV/2-l+i/)! . . 

My) = W 2-l)!^ o {N-l + 2v)\ HN - l + 2v{y) 



(2- 



7r l/2 e -a 2 /2 - 



{N/2 - 1)! ^ 2^+ 2 '(iV/2 + 0! 



1 



H N +2i{y) 



(in these formulas it is assumed that N is even). The task now is to show that with 
J = (s, oo) and the soft edge scale s h-> \/2iV + s/2 1/2 7V 1/6 , the large iV limit of 
(TO!) with (0(2) given by ([121) is e q ual to We know from [Sos02j (see |BF03] 

for a restatement of this) that for this task it is sufficient to show 



lim / dxi 

A->oo JV2N+s/2 1 / 2 N 1 /e 



dx o odd{GOBN)2 (xi x ) 

2N+s/2 1 / 2 N 1 / 6 v 1 



(2.9) 



dx\ 



, odd(OEsoft) 2 , \ 



The starting point of the derivation of (|2.5|) given in [ForOOj is very similar. 
In the notation of (|1.15|) . the analogue of (12. 7|) is 



(2.10) 



^odd(LOE„) ^ ...,£„) = det K N ^!(xj,Xk) + (.92(2^)) 1/2 L a N _ 1 (x j ) 

(N-iy. 



' 2*-2((JV - 2)/2)!(a/2 + (TV - 2)/2)! 



B^Xk) + B 2 {x k ) 



j,k— l,...,n 
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Here Kn-i(x,ij) refers to (|1.6|) with gi{x) = x a e x and (assuming N is even) 



!/=(iV-2)/2 v 



D , , ^ 2 2 >/2 + ^)W . 1/2 

(2Q! , 
2 2l l\(a/2 + l)\ 



r9in ^ ^i (( a -l)/2)! 2 (a/2)! 2 (2Q! 1/2 

(2.11) £? 2 (y) = 2 — ^ 920ir»/9^/V (g2(y)) L2i(y) - 



l=N/2 



Here the task is to show that with J = (0, s), and the hard edge scale s >— > s/AN, the 
large AT limit of (|1.4[) with /j( 2 ) given by (|2.1[) is equal to (|2.5p . For this, according 
to [Sos02 , it is sufficient to show 



s/4JV m/4/V 2 

lim / dxf- dx n p°^ d(LOEiv) (xj., . . . , x n ) 



(2.12) =y dXf-J dx n ^ddCOEhard) ^ . ..,£„). 

A mechanism for the validity of (|2.9|) and (|2.12|) is the uniform estimates 
(^I7i) "<) d(LOEN)2 (V2iV + Xl /2^N^, ...,V2N + x n /2^N^) 

(2.13) = pffi OEmft)2 (xi, ...,x n ) + o(l)R s n ( Xl , ...,x n ) 
and 

{^)y^ W \*rm...,x n /W) 

(2.14) = p^OEh-rd)" _ + o(1)i? h (a;i) ^ Xn) 

where o(l) refers to the dependence on N, while R n is integrable on the domain in 
question. In |ForOO] only the leading term in (|2.13[) and (|2. 14[) was computed, so 
in particular (12. 3|) and (|2.5|) were not rigorously established. This "technical issue" , 
essentially asking for uniform estimates of the infinite sums (|2.8|) and (|2 . 1 [> . is the 
reason that [ForOOj was not submitted for publication. 

2.2. Special superimposed ensembles. With the task of providing uniform 
asymptotics of (|2.8[) and (|2.1[) being technically difficult, necessity dictates seeking 
an alternative strategy. For this one should bring to the fore the notion of univer- 
sality, which tells us (for example) that there is nothing canonical about the finite 
N ensemble (|1.14|) in regard to studying the limiting distribution odd(OEsoft) 2 . 
Thus for a general weight function gi(x) 

2n 

E l (n; (s, oo); odd(OE JV ( 5l (a;)) U OE N ( 9l (x))) = ^ E^2l - 1; (s, oo); 9l (x)) 

1=0 

(2.15) x (E^l; (s, oo); 9l (x)) + E x (l - 1; (s, oo); fll (x))) 

which follows immediately from the definition of the superimposed ensemble. Hence, 
the universality of the gap probability in the superimposed ensemble is a conse- 
quence of the universality at the edge of the ensemble OEjv(<7i(x)), which is known 
from IDG06] . In particular, instead of considering the soft edge scaling of the 



GAP PROBABILITIES 



9 



superposition of Gaussian orthogonal ensembles, we may just as well consider the 
soft edge scaling of superimposed Laguerre orthogonal ensembles 

(2.16) odd^E^O^e-^ 2 ) U OE^t/e" 1 / 2 )) , 

and more particularly this for any one value of a. So the question now is, amongst 
the ensembles (|2.16[) is there a value of the parameter a for which the correlations 
have an explicit form which is easier to analyze that that in (12. 7j) ? 

That there is a special ensemble amongst (|2.16j) is seen from a theorem in 
[FROlj . This theorem classifies all continuous weight functions g\ , 52 such that 

(2.17) evenf OE„( 5l ) U OE n ( 9l j) = UE„( 32 ) 



where UE„(g2) refers to Ijl.lj) with (3 = 2 (here the U denotes unitary symmetry 
and even refers to the labelling of the eigenvalues countered from the right, which 
is the soft edge) . In fact up to a fractional linear transformation there are only two 
weight functions with this property, 

(2.18) I {{1 _ x) (a-W {l _ xn _ 1<x<1 . 

Hence of the superimposed Laguerre ensembles in (|2.16p . the case a = is dis- 
tinguished by the property (|2.17j) . In keeping with this special feature, for the 
correlations of every odd labelled eigenvalue as required by (|2.16|) . it allows the 
structured formula [FR04] 



(2.19) p°^ {WE - )2 ( Xl , . . . , x n ) = det [ - A f K K N (x„u) du 



Ox 



where refers to (II. 6|) with g^x) = e~ x , Pj(x) = L®(x) (here L®(x) denotes the 
Laguerre polynomials with parameter a — 0). 

One sees immediately that the structure exhibited in (]2 . 19[) gives a much cleaner 
expression than that exhibited by (|2.7|) . The task is to compute uniform asymp- 
totics of this under soft edge scaling, which for the a = Laguerre ensemble is 
obtained by replacing coor dinates [For93j x h-> AN + 2(2iV) 1 / 3 X then taking the 
limit N — > 00. It is a straightforward exercise using the uniform estimate |OIv74j 

r_iyv 

(2.20) e-*'*L° N (x) = 7^jb Ai W + O^WN- 1 ' 3 ) 

where x = 47V + 2(27V) 1 / 3 £, valid for t € [to, 00) to obtain the uniform asymptotic 

expansion 

(2.21) 

2(2iV) 1 / 3 J fi: A r(4iV + (2A0 1/3 s, 4iV + {2N) 1 'H) = K soh (s, t) + 0(e- (s+t) )0(A^ 1/3 ) 

valid for t, s G [to, 00). Further, this expansion can be differentiated with respect to 
t or s. However, as written in (|2.19[) the argument u in Kjsf(xj,'u) takes on values 
which are to leading order in [0, 4JV] instead of [4iV, 00) as in (|2.2ip . As noted in 
FR04] Lemma 13], this can be circumvented by using the identity 

px poo J poo 

K N (y,u)dy=(-l) N - 1 e - u / 2 —L° N (u)du- K N (y,u)du. 
,/n ./„ du .L, 
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Use of (|2.20j) and (12.2ip then gives the uniform asymptotic expansion 

r4N+2{2N) 1/a X 

2(2N) 1/3 / K N (4N + 2(2N) 1/3 Y,u) du 



(2.22) = / Ai(«)du- / K sott (Y,u) du + 0(e- {x+Yy )0(N- 1/3 ), 

Jy Jx 

which furthermore remains valid upon differentiating with respect to X or Y . Not- 
ing from the integral form of the kernel (jl.9[) . 

/>oo 

K sott (x,y) = Ai(x + u)Ai{y + u) du 
Jo 

that 

f K soft (Y,u)du = K solt (X,Y) + AUY) f M(t)dt, 
and substituting this in (|2.21[) then substituting the result in (|2.19p we deduce 
(2(2N) 1 / 3 ) n p° ( ^ (LOE ° N) \iN + 2(2iV) 1 / 3 a; 1 , ...,4N+ 2{2N) 1 ' 3 x N ) 

(2.23) = p°^ {OEsott)2 (x u ...,x n ) + e"^+-+^0(N-^ 3 ) 

where p° dd ( OEsoft ) j s specified by (|2.4|) . ft is immediate from this that the analogue 
of (12. 9p holds true. But according to |Sos0 2] the latter is sufficient for the validity 
of the limit formula 

(2.24) lim £° dd ( LOE «) 2 ((4/v + 2{2N) 1 ' 3 s, oo); £) = ^(OEsoft) 2 ^ ^ ^ 

N—*co 

with £: odd (OEsoft) 2 ^ oo); £) specified by (|2.3|) . Hence, for the soft edge scaling, our 
sought identity has been established, providing us with a rigorous justification of 
the identity 1(1712"]) . 

It remains to establish a similar limit formula for the hard edge. Of the two 
special pairs of weights (|2.18p , the second pair near x = 1 exhibits a general hard 
edge singularity. Recalling that the ensemble OE N ((l—x) a (l+x) b ) with — 1 < x < 1 
is referred to as the Jacobi orthogonal ensemble JOE^ 6 , the left hand side of (|2.17p 
for the second pair in (|2. 18[) refers to the even labelled eigenvalues counted from the 
right in the ensemble JOE^~ 1)/2 '° U JOE^~ 1)/2 '°. Analogous to ([2721]) the n-point 
correlation function for the odd labelled eigenvalues of this ensemble is given by 
the structured formula [FR04, eq. (2.16)] 
(2.25) 

odd(JOE N ) ( Xl ,...,x n ) =det -- — (1-Xj) (l-Xj)K N (xj,u)du 

y ' L OXj J_i lj,k=l,...,n 

Here K N (x,y) = K N (x,y)/((l - x)(l - y)) 1/2 where, with p( a ^(x) denoting the 
Jacobi polynomials, Kff(x,y) refers to (|1.6p with g-2{x) — (1 — x) a and Pj(x) = 

Pf'°\x). 

In (|2.25p the integration variable u is not confined to the neighbourhood of the 
hard edge. To avoid this potential problem, we make use of the identity [FR041 
eq. (3.69)] 

(2.26) / K N (x,u)du = -2 [ K N (-l,u)du 
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where we regard Kn(— l,it) specified by the final form in (11. 6|) . supplemented by 
the special value Pj a '°\-1) = (-T) j . 

The hard edge scaling limit in the neighbourhood of x = 1 requires replacing 
the coordinates x i— > 1 — X/2N 2 then taking — > oo. To analyze the latter limit, 
we make use of the uniform asymptotic expansion |Sze75j 

(sm|)°(cos|) P^(cos0) 



(2.27) = n~ a v , V — J a ((n + (a + 6 + l)/2)0) + 6? a+2 0(n a ) 

m V smw 

valid for < 9 < c/n, (c > 0), which furthermore remains valid upon differentiation 

with respect to 6. With use made too of (|2.26[1 . this gives 

(2.28) 

/y p°° / 1 \ 

K N (x,u)du ^_ X/2N2 = -X 1 ' 2 / v-^ 2 K^ d (X,v)dv+0(-)0(l), 

which furthermore remains valid upon differentiation, with the dependance on X, Y 
again O(l). Using the integral form of the kernel (|1.16|) . 



If 1 _ 

K h ^ d (x,y) = - J a (Vxl)J a (Vyi)dt 
4 Jo 



we have that 

I"" v ~^K h3lA {X, v) dv = K h3xA {X, y) + J *(^ r J a (t)dt. 
dX J Y 2yY J^/y 

Substituting this in (|2.28|) , substituting the result in (|2.25[) , and recalling (|2.6[) we 
conclude 

(2.29) =^< OEh " d >"( :Cl x n ;(a-l)/2) + 0(l)0(l) 

where p° dd ( OEhard ) j s specified by (|2.6j) . From this uniform asymptotic expansion 
it is immediate that the analogue of (|2.12p holds true, and again by appealing to 
|Sos02] the latter is sufficient for the validity of the limit formula 
(2.30) 

lim ^ dd ( JOE «" 1,/2 '°) 2 ((! _ S /2N 2 , 1); = £ odd ( OEhard ) 2 ((0, s); £ (a - l)/2), 



N- 



with £;odd(OEhard) 2 spe cified by lj2~5)) . This is the sought limit formula for the hard 
edge scaling, complementing (|2.24l) for the soft edge scaling, and providing us with 
a rigorous justification of p. 131) . 
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